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Abstract 


The procedure developed to extract average source direction and average 
source size from spin-modulated radio astronomy data measured on the IMP-6 
spacecraft is described. The method is quite general and fast and is easily 
implemented. Because all measurements are used, rather than Just finding maxima 
or minima in the data, the method is very sensitive, even in the presence of 
large amounts of noise. The technique is applicable to all experiments with 
directivity characteristics. It is suitable for onboard processing on 
satellites to reduce the data flow to Earth. The application to spin-modulated 
non-polarlzed radio astronomy data is made ia uetall and includes the effects of 
noise, background, and second source interference. It is shown, for example, 
that measurements made using low-resolution dipole antennas with a sampling rate 
of about 2 measurements per spin period yield a determination of average source 
position with 1 deg accuracy for signals greater than several times background. 
This can be achieved utilizing only 25 samples from a receiver such as the 
Goddard Space Flight Center Radio Astronomy Experiment on lMP-6. For large 
sources such as solar bursts, the average source size can also be determined to 
a significant accuracy. In addition, the method permits the accurate restora- 
tion of the signal by removing the modulation. The analysis has been tested 
with computer simulated data and the results agree with analytic predictions. 
Applications of this method with IMP-6 radio data have led to a) determination 
of source positions of traveling solar radio bursts at large distances from the 
Sun, b) mapping of magnetospheric radio emissions by radio triangulation, and 
c) detection of low frequency radio emissions from Jupiter and Saturn. 


1. INTRODUCTION 


Many experimental measurements deal with physical quantities having a 
magnitude and direction (vector), such au fields (electric or magnetic) and 
fluxes (particle, photon). Most experiments are designed with a directivity 
(beamwidth) which is necessary fcr the determination of source direction and 
angular variation. Groundbased experiments emphasize narrow directivity and 
high sampling rate during an angular scan in order to resolve and measure 
sources of angular size comparable to the experiment beamwidth. However, for 
spacecraft measurements, the design of an experiment is greatly restricted by 
weight and power limitations. A narrow beamwidth for a given aperture size 
often implies reduced sensitivity; a high sampling rate to allow several samples 
per beamwidth requires a high telemetry bit rate. As a general rule, it is 
necessary that an experimenter not restrict beamwidth more than is required for 
the phenomena measured in order to conserve both experiment weight and telemetry 
resources. 

In addition, there may be other factors affecting the measurement process 
which are even less under the experimenter's control. For example, the spin 
rate which determines how fast sources are scanned is a compromise between the 
needs of all experimenters and the requirements of spacecraft dynamics. In view 
of the conflicting requirements between spacecraft experiments and spacecraft 
operations, it becomes increasingly important to extract from the spacecraft 
data in an efficient way whatever information is present. 

In this paper, we describe a data analysis technique developed for and first 
used with IMF-6 radio data (Fainberg et al,1972) which has applicability to a 
wide variety of spin modulated measurements. The technique is very fast, so 
that it is possible to process a large amount of data. It has great flexibilty 
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and sensitivity. In addition, with the source position and size determined, it 
is possible to accurately remove the effects of spin modulation and thereby 
determine the nature of the data in the absence of spin. 

In Section 2, the analysis technique is discussed for general spin modulated 
experiments. In Section 3, the application of the analysis to spin modulated 
radio flux measurements is made in detail. Formulas are developed to account 
for effects of background, noise, and other sources on the measured source 
position and size. In Section 4 these techniques are applied to simulated radio 
bursts, with results agreeing quite well with analytic predictions. In 
addition, we show some examples of actual data processed frcq the IMP-6 
spacecraft. Section 5 is concerned with future applications and extensions. 

The purpose of this report is to discuss the data processing technique 
developed for the GSFC radio astronomy experiment on IMP-6, and to illustrate 
the technique with some applications to IMP-6 data. For an excellent review of 
low frequency satelite radio astronomy techniques in space, particularly with 
respect to questions of polarization , see Lecacheux et al [1979]. For a 
general review of spacecraft observations see also Steinberg [1979]. 

2. GENERAL METHOD 

Initially we consider the case of a distant point source located in the spin 
plane of the experiment. Later we relax this restriction and show how an 
extended source out of the spin plane can be treated. 

The spacecraft experiment has a normalized directive response or gain G((|i) 
which is a function of direction (j) in the spin plane. For convenience we 
assume the maximum gain at (])= 0, viz G(0)= 1. The gain function is determined 
by the design and calibration of the experiment. In all experiments, since G(i<i) 


is periodic in ^ (modulo 2ir), it can be expressed as a Fourier series, 

G( 4 >)= z (bj^cosk 4 i + ajjSink^) (1) 

where aj^ and bj^ are the Fourier coefficients obtained from G(4i) in the usual 
fashion. Very often symmetry considerations reduce the number of Fourier 
coefficients. For example, if the gain is symmetric about the direction of 
maximum gain, i.e. if G(<i)= G(-4>), then aj^= 0 for all k, and 

G((ji)i z bj^cosk^i (2) 

An example of a specific gain function is that of a dipole antenna with 
length short compared to radiation wavelength: 

G(i}i)= |(1+cos24>) (3) 

For antenna lengths comparable to wavelength, additional Fourier terms become 

significant. This case is treated in detail in Section 3. 

For the general case let a distant small source of strength A be located at 
4 . Then the experiment response at time t will be 

5 

R(t)= G(4.g-a)t+u)tQ)A (4) 

where 

( 0 = spin angular velocity= 2ir/spin period 

tps time when maximum gain is oriented in the reference direction (<((=0) 

The source may be time varying in amplitude and position so that 

R(t)= G(<^g(t)-ait+(i)tQ)A(t) (5) 

where A(t) is the experiment response to the source at unity gain in the absence 
of spin. Since measureraents are made at discrete times t^, tj^_^^,..., the 
observed quantities form a time series R(tj^), R(t^_^^),... The general problem 
we consider is the following: given this series of observable quantities (R(t^)} 
and the gain function of the experiment G(4>), how can we determine the source 
position and the source strength as functions of time, and what is the precision 
of the method taking into account noise, background, and the presence of other 


aouroes? 


Tha procedure we use ia baaed on Fourier techniques in which we choose a 
processing window of N consecutive observations and form the sums and over 
N points where 

Z R(tj_)oos(wt^-wtQ) 

S^s S R(tj. )sln(tot^-wbQ) (6) 

It is shown in Appendix A that these quantities reduce to 
1 1 

I b^cos<i)g+5 a^simtig + small terras 

1 1 

S^= a^cositig+j b^simi)g small terms (7) 

Therefore 


tani|>„= 

s 


b^S^+a^C^ 

b^C^-EiSi 


( 8 ) 


The resulting source angle <tig refers to the time of the midpoint of the 
processing window. The source angle at one sample time later can be obtained by 
shifting the times in Equations 6 by one sample time. (Note that this is easily 
done by adding one new term and subtracting one old term from the sums in 
Equations 6) . In this way the evolution of source position with time can 
easily be calculated from experimental data. 

Note that if the gain function (Equation 1) has significant Fourier terms of 
order k, then we also have 


tank(j)„= 




( 9 ) 


where 

Cj^= E R(tj^) COSk(o)b^-utQ) 

Sj^= Z R(t^) sink(wtj^-uitQ) (10) 
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This allows an independent determination of source position utilizing each 
set of Fourier gain coefficients. 

3. APPLICATION TO RADIO ASTRONOMY MEASUREMENTS 

In this section we consider in detail the processing technique described in 
Section 2 applied to non-polarized radio data measured with a dipole antenna. 

In Section 3a the dipole gain function is discussed and the nature of the signal 
observed from a distant point source is introduced. In Section 3b we consider 
the effect of several point sources, discuss the response of point sources out 
of the spin plane, and then consider the case of a distributed source out of the 
spin plane. For these sources we introduce a quantity called the modulation 
index, specifying thi amount of observed modulation. In Section 3o we give the 
procedure for determining the modulation index' and source angle from the data. 
Section 3d is concerned with the accuracy of the source angle and modulation 
determination taking into account signal to noise ratio, effects of background, 
and the size of the processing window. In Section 3e we discuss the procedure 
for removing a modulated background from sporadic events. 

3a. P ix>.o l e P ain Funatign 

We consider the normalized gain function of a dipole antenna which spins 
about an axis perpendicular to the dipole elements. The angle ij> is measured in 
the spin plane which also contains the dipole. (The case of a dipole nob 
located in the spin plane will be treated in a subsequent paper) . Angles are 
measured from the direction of maximum gain, which is the direction in the spin 
plane normal to the antenna elements. From Equation 2 we have 
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G(4>)s S bj^OQsk.^ (11) 

since G((i>)s G(-4i) , This antenna possesses another symmetry, G(t)= G(4i+180), 
which implies bj^s 0 for k= 1,3,S,.*> Iherefore for a dipole the most general 
gain expansion is 

G(i|>)s bQ+b2OOs2$+bi|0Osl|i^+, . . (12) 

For antennas short compared to the wavelength of the radiation measured, only 
the first two terms of this series are significant and 

G((^)= g + gcos2<t> (13) 

An expression for the gain of a dipole antenna which is not restricted to very 
short lengths is given by Jackson [1962]: 

G(4>)s Gq[ cos|(kdsin^)-cos|kd ]^/ cos^iji (14) 

where 

k= 2ir/X 

X= wavelength of radiation 
ds total length of antenna 
Gq= [1-cos|kd]’^ 

Note that this reduces to Equation 3 for kd << 1. 

We can expand this general dipole gain function in the form of Equation 12 
by standard B’ourier techniques. We have done this by numerical integration and 
have plotted the coefficients b^, bg, bj^, and bg in Figure 1 as a function of 
normalized antenna length. It is clear that Equation 13 is an accurate 
representation for the gain for wavelengths larger than half the dipole length. 
For shorter wavelengths the terms bg, etc may be large enough to be useful 
in the procedure given in Equation 9. 
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3b> Nature of Obaerved Signala 


To illustrate the development of the techniques in detail, wo consider the 
range of situations where Equation 3 is sufficient. In most of the development 
it is possible to generalize the treatment to include the higher Fourier terms 
and to include corrections for bg/bg 4 1, 

For the short dipole the observed signal from a distant point source in the 
spin plane is given by Equations 5 and 3 , 

R(t)s 2 Cl+cos 2 (uMg)]A(t) (15) 

where A(t) is the signal observed by a unity gain antenna in the absence of 
spin. For most mses of interest A(t) and t^Ct) are slowly varying functions 
compared to cos2i»it. These cases are illustrated in examples in Section 4 where 
the spin modulation is rapid compared to the slower variations in actual data. 
When A(t) has more rapid time structure the effect will be seen to be somewhat 
equivalent to a poorer signal to noise ratio in the determination of ♦ using the 
techniques proposed. We write Equation 15 in the form 

R(t)= [1+ cos2(ut-4'^)]A^ (16) 

keeping in mind that source strength A^ (which now includes a factor of ^) and 
source direction may be slowly varying in times. 

In the remainder of this section we will show that any distribution of 
sources results in an observed signal of the form 

R(t)= [1+ a cos2(u)M)]A (17) 

where A is a measure of the amplitude of the source, t is the arrival direction 
of the midpoint of the source projected onto the spin plane, and a is defined as 
the modulation index such that 0<a<1. To show this we consider 3 cases: two 
point sources located in the plane of spin; one point source located out of the 
spin plane; and a uniformly distributed conical source. 




For two point souraes looatsd at and <>g we have fon the observed, signal 
R^^Ct)* R^(t)+Rg(t) (18) 

where R.j(fc) and RgCt) are given by Equation 16, By standard trigonometric 
procedure it is possible to write this in the form 

Rj(t)» Cl+a|,oos2(wt-^j)3Aj (19) 

where is a source angle representing a weighted moan of (j>.j and (^g, and is 
the sum of the Individual source strenths A.| end Ag. 

Here «^<1 and is a positive quantity which we define as the modulation 
index., When the two point sources are in the same (or opposite) direction, a^sl. 
When two equal sources are at 90**, c^.sO (no modulation). The modulation index 
is than a measure of the separation of th two sources. Rather than giving 
detailed expressions for and for two point sources now we defer this until 
later in this section when we consider the sum of two point sources not 
restricted to the spin plane. 

Case 2 

We consider a point source located out of the spin plane. It is clear that 
we cannot have c= 1 because the otaerved intensity R(t) can never be zero since 
the null of the Jintenna no longer passes through the entire source. This case 
is treated exactly in Appendix B and we again can express the result in the form 
of Equation 17 where <}! is the azimuthal angle of the source in the spin plane 
and the moduation index is given by 

(t- sin^0/(2-sin^S) (20) 

where e Is the angle of the point source to the spin axis. We note that here 
again the modulation index is between 0 and 1. However, lack of full modulation 
in this ease Is due to a point source out of the spin plane. 

The modulation index a is a quantity which we can determine along with source 
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direotion from the obaerved data. However, it is clear that a unique oonolusion 
from a is not possib)*,* laok of oompXete modulation oan be due either to a point 
source out of the spin plane or to distributed sources. Ihe range of ambiguity 
is restricted, however, and the modulation index does contain important 
Information, We can explore this aspect by treating next a very general oasei 
that of a distributed uniform source centered out of the spin plane. 

JSaflgJi, 

We consider a uniform conical source of angular width y, located with its 
axis inclined to the spin axis at an arbitrary angle 0. The signal observed 
from this source has been expressed in closed form in Appendix B, We find that 
R^b [1+acos2(ut-i^)]A (21) 

Where the modulation index is given by 

as sin^0/(K-sin^0) (22) 

where 

K= [84-2oosy( Ucosy) 3/(3oosy( 1+oosy) 3 
Y= half total angular width of cone 

We note that for 0s 90" and for ys 0 this expression describes the case of a 
point source in the spin plane (os 1). For ys o this case refers to a point 
source out of the spin plane, and Equation 22 reduces to Equation 20. 

For a given a, only a certain range of source widths and source inclinations 
is possible. Figure 2 presents a graphic description of this range for various 
values of a derived from Equation 22. In this figure the vertical axis 
represents the elevation of the source from the spin plane. The horizontal 
scale is half the total width of the uniform conical source. Once a value of 
modulation index is found the allowed range of source widths and source 
inclinations is given by the appropriate curve in Figure 2. The closer the 
modulation index is to unity the less the range of ambiguity. 
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3o 


We sbsrt wiWi the tBoab gineral form of obaerved signal from a short dipole 
antenna located In the spin plansi ineludlng now an unmodulated background B, 
RCt)* i’l4aco82C«t-t)3A + B (23) 

where R(t) is observed only at discrete times tj^, We wish to determine from 
{R(t^)) alone the quantities n, <|i, and A which, in general , are slowly varying 
in time eompared to the modulation imposed by the spin. We assume B can bo 
determined at a time when the source is not present , Following the procedure of 
Section 2 we fom the sums 

Cgs R ( b|^) eoaSfcb 

SgS S RCtj)ain2«tj_ (2R) 

and also 

QqS 1 R(fcj) 

where ws ‘dn/apin period . 

In order bo explore the moaning of Sg, Gg, and Cq, we will assume Equation 
23 and insert this in Equations in order to find expressions for a, and A 

in terms of the observed data iR(tj^)} only. In addition, it is important to 
note that A and B are noise quantities that are sampled by an instrument with 
certain physical parameters. It can be shown (Kraus, 1966) that the rms 
fluotuation AA renultlng irom a radiometer measuring a noise signal of average 
power A is 

AA« A//(0t) (25) 

where 6 is the predoteetlon bandwidth of the radiometer and t is the post- 
detection Iniegrabior time, fJimilarly 

b/V(6t) (26) 

where AB io the rmo fluctuation in the roeasurod background power B. Using 
Equations 13 and 1^ we find (Appendix A) that 
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SgS gNaAsina^ iO(AS2) 

Cgs |NaAoos2* ±0(AC^) (27) 

Cq= NA+NB ±0(ACq) 

where 0(ASg)i O(ASg) » and O(ACq) are anall terms of the order of the standard 
deviations. Neglecting the small terms, Equations 27 can be solved for a, i)>, 
and A: 


tan24>s S 2 /C 2 

(28) 

as 2/[S2+C2 ]/[Cq-NB] 

(29) 


or finally 

Z R(t j^)sin2i>)t^ 
tan2(|)- j, R(t^)cos2(i)t^ 

a= 2/C( £ R(t^)sin2a)tj^)^+(r R(t^)sln2ut^)^]/£ [R(t^)-B] (30) 

As CR(t^)-B] 

which expresses the desired quantities in terms of the observables R(t^) and the 
spin period. We note that Equation 30 is a special case of the general solution 
given in Equations 9. 

To reconstruct the signal that would be observed by an antenna with no spin 
there are two possible approaches. One method is to use Equations 30 and 
multiply by (1+a) in order to insure the appropriate normalization 

A(tj^/2^= r(1+c«) 2 [R(t^)-B] (31) 

This is quite useful but it does involve a smoothing over the window length of N 
samples. Another method is to use Equation 17 which yields 

Ao(ti)= (1+o)[R(t^)-B]/[1+ocos2(ut^-4>)] (32) 

where a and <b are determined in Equations 30. This function, which we call the 
despun function, does not involve a smoothing of A(t) over the window (except 
for the determination of o and (|i) and more accurately follows variations with 
time, but at the expense of more noise. In the next section we consider how the 
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accuraoiea of 4', and A are affected by the error terms in Equations 27, 

3d. AsguraoY.. . of. ...Method ; Unmod ul ated Background 

We consider first the error in ^ due to the terras ASg and ACg. If AA and aB 
are the rms fluctuations of signal and background respectively, and A 41 and Aa 
are the rms values of errors in 4 and a respectively, we have shown in Appendix 
A 

A$g= V[1+^a^+(1+a)^ AA^/AB^]aA/[A/( 2N)] radians 
Using Equations 25 and 26 we get 

A4»g= “/[1+^a^+(1+a)^ A^/B^3 /[/(2Net) ] radians (33) 

For large A/B 

Aa= /C2-.3a^/2+c»'*/23//(NpT) (34) 

Several implications of these results are clear. Once the signal to background 
ratio A/B is greater than about 6, the background has a small effect on the 
accuracy of 4'. On the other hand, if the signal A is long lasting, one can 
reduce the errors in 4> and a by increasing N, which is effective even for very 
poor signal to noise ratios. 

Consider the case where A/B s 5 and 1//(et) = O.O 9 (value of Pt for the GSFC 
radiometer on lMP-6) . Then for N= 25, A4i“ 1 deg and 6a« 0.02. This indicates 
the great sensitivity of this data processing method. 

There are several considerations governing the choice of N. Large processing 
windows decrease errors by 1//N, However, they also smooth out structure within 
the time period of N samples. This effect is very similar to the use of a 
rejection filter to get rid of noise signals, but at the expense of losing 
signal components within the same frequency range. Another consideration is 
that the number of samples N be such that an even sampling of phases occurs in 


the sums for appropriate oanoellatlon. We have achieved good results with N as 
small as 12. Very often there is an aliasing between the sampling time and the 
spin rate. This occurs when the sampling time is such that the detector has 
rotated a nearly integral multiple of half spins to an equivalent position only 
slightly different from that of the preceding sample. In such oases, it is 
useful to have N cover an integral number of aliasing periods. We have also 
found it useful to do an additional smoothing ( running average) on the derived 
values of ^ and a to reduce any remaining residual modulation caused by a 
non-uniform sampling of phases. 

3c Jm. Source Effects 

We consider the effect of having two sources present with arbitrary 
strengths, positions, and modulation indices. Let the stronger source be 
located at =0 and the weaker source Ag be located at Let represent 

the observed sum of the two sources. We have 
R^(t)= A^[1+a^cos2wt] 

Rg(t)= A2[1+02°o®2(a)t-(ti2i)] (35) 

R^(t)= A,p[1+a,pC032(a)t-(|ij)] 

and 

Rj(t)= R^(t)+R2(t) (36) 

Then 

A,pZ A^+Ag 

a,p= a^/[1+2cos2(|»g^ Ag02/(A^a^)+ (A20g)^/( A^a^)^ ]A^/(A^+Ag) (37) 

tan2i|)^s (A2a2sin2(j)2^)/(A^o^+A202sin2i|»2^) (38) 

Note that 4>.p, which is the measured shift ft?om source 1 position due to the 
presence of the weaker source 2, is zero when = 90 “, i.e., when the two 
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sour’oes are at right angles. When «2 = ° source 2 is isotropic, o® and 

a |/( l+Ag/Ai ) ' (38) 

The behavior of Equation 37 foi^ is plotted in panel A of Figure 3 for a 
range of source strengths as a function of source separation. Panel B shows the 
fractional change of modulation index under similar conditions. It is clear 
that a weak source does not appreciably affect the der ived position or 
modulation index of a strong source. 

3f .Removal of Modulated Background 

In most cases, transient events such as solar bursts or planetary emissions 
are observed simultaneously with a modulated background component. If the 
background characteristics ( amplitude, direction, modulation index ) are 
relatively constant, then they can be determined before and after the transient 
event and corrections can be made to obtain the burst characteristics alone. 

Let Aq, Oq, and i})q be the average intensity, modulation index, and position 
of the background and let A.j , and (f.j he the measured average intensity, 
modulation index, and postion of the transient event plus the background. We 
wish to find A^, a^, and ij)^ characterising the source alone. These quantities 
are related by 

Ag[1+a^cos2(a)t-it)g)]s A.j [ 1 +a.|C0s2(wt-<j) .j )-AqC UaQCOs2(tot-(t>Q) ] (40) 

After simplifying we obtain 

Og= /[ (a.jA.j)^-(aQAQ)^-2a.jA^aQAQCOS2(<f^-'t>o)]/(A^-"AQ) (4l) 

tan2(j)g= Ca^A.jSin 2 ()).j-aQAQSin 2 (l)Q]/[a.jA.,cos 2 ()).j-aQAQCOs 24 )Q] 

By Standard propagation of error techniques it is possible to derive expressions 
for the standard deviations of these source determinations in terms of the 


standard deviations of the background (0) and the background plus signal (1) 
parameters. However, the expressions are long and the calculations are best 
performed during the data processing. 

4. APPLICATION TO SIMULATED DATA AND TO REAL DATA 

In order to verify the accuracy and features of this processing method, we 
have found it very useful to generate data with characteristics similar to 
actual data. For the simulated data we specify the amplitude, position, and 
modulation index and then see how well the processing method works. 

In Section 4a we consider the technique applied to a signal with an abrupt 
start in the presence of background noise. We show that it io possible to 
determine the direction of arrival of signals having an inter tjity only 0,1 that 
of background. In Section 4b we consider the technique applied to a simulated 
solar burst. In Sections 4c, 4d, and 4e illustrations are given of spacecraft 
measurements of real data from solar bursts, Earth emissions, and other 
planetary emissions. 

4a. Simulated Step Signals 

In this section we illustrate the results of the processing technique 
applied to examples of simulated data consisting of a step signal in the 
presence of background. The input data are computer generated with noise 
fluctuations having similar statistics to those of data obtained with the 
Goddard Space Flight Center (GSFC) radio astronomy experiment on IMP-6 (see 
Brown [1973] for an experiment description). For all examples Af/f= 0.09 where 
Af is the standard deviation of simulated receiver power with a Gaussian 
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disfcribufclon. A data sampling rate of 11,7 samples per minute and an antenna 
spin rate of 5*39 rotations per minute are used. These values are similar to 
those in the lMP-6 experiment and result in observed beat or aliasing periods 
of 1.05 rain (12.3 samples). However, the data-prooessing technique is not 
dependent on any particular sampling or spin rate as long as a relatively 
uniform sampling of phases is used in the processing window. For the oases 
treated here, the best processing windows contain the closest number of data 
points to integral multiples of 12.3. 

In Figure A panel A is a plot of the input data (with no spin) as a function 
of time. Before t=0, just the background is present. At t=0, a step signal is 
switched on which Is 10 dB above (10 times) background. Panel B of Figure 4 
shows these input data with modulation caused by the spacecraft spin. The 
signal characteristics for this case are signal amplitude =10 dB, modulation 
indexs 0.75, and direction of arrival = 0°. The background characteristics are 
amplitudes 0 dB and modulation index = 0 (non-directional) . Using each 1 minute 
(12 samples) of signal in panel B as input data processing technique 

yields the results plotted in panels C, D, and E. Panel E shows the derived 
direction of arrival (with a 12 point smoothing) as a function of time. The 
curve marked "total" is the result of Equation 30 applied to the data in panel 
B. Prior to time t = 0 min, the derived positions exhibit a random walk as 
expected for an unmodulated signal. At t = 0 min the computed position shifts 
towards 0°, which is the true position of the input signal. In panel D, the 
curve labeled "total" results from Equation 30 applied to the input data. 

During a ten minute period prior to t = 0, the characteristics (mean intensity, 
position, and modulation index) were determined for the background using 
Equations 30; the values obtained were then used in Equations 4l to subtract the 
background effects. It is seen that for this case a significant correction 
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towards the true value, a x 0.75» occurs , indicating a proper subtraction of 
background, Panel C is a {d.ot of the derived mmodulated input signal using the 

values of position and modulation index plotted in panels D and E. Upon 

comparing panel C to the actual input in panel A, it is seen that the signal 
restoration is very good, except near tx 0, which is at the region of 
discontinuous change. 

In Figure 5, a case is exhibited where the signal is 20 dB above (100 times) 

background. It is clear that the signal is so large that there is little need 

for background removal. Figure 6 represents a case with smaller signal 
intensity now equal to background. Even in this case the true position of 0° is 
well measured. The modulation index exhibits more scatter because of the lower 
signal to noise ratio, but, nevertheless, the derived modulation index is an. 
approximation to the input value of 0.75. The errors are similar to those 
predicted by Equations 34 and and error analysis of Equations 40. 

The above simulations were run with a processing window N of 12 data points. 
To illustrate the power of this signal processing technique for longer 
processing windows, a value of ns 123 points was used in the case shown in 
Figure 7> In this example the signal amplitude is again equal to the background 
as in Figure 6, but now is processed with a 10.5 min window. It is evident that 
a greatly improved determination of the true source position and source 
modulation index are obtained with this longer window. It is then possible to 
determine characteristics of even smaller signals. Figure 8 shows a case for 
the signal intensity only 0.1 that of the background. This input signal is 
not visible at all in panels A and B. However, the method yields a respectable 
determination of position (0* ) in panel E and modulation index (0.75) in panel 
D, even for this case where the signal is burled in the noise. 

As a final example we show a case where the background is directive and 
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therefore modulated before and during the time of the signal. For this case, 
shown in »;’igure 9, the baol^round charaoteristios are modulation index * 0.25 
and direction of arrival = -60". The step signal starting at t= 0 has a 
modulation index = 0.75, direction = 0", and signal amplitude equal to the 
background amplitude. It is clear from panels E and D that with background 
removal, the processing technique yields a very accurate determination of the 
true signal direction and modulation index. 

^jBiulat..ed..&Qlar Burst 

A very common form of sporadic radio event measured by spacecraft 
experiments is the rapid drift (type III) solar radio burst. We found we can 
sJ.mulate such a burst by a simple function of the form 

A(t)=t“ exp(-at) (42) 

This function has a maximum (n/a)*^ exp(-n) at time = n/a. If the maximum 

amplitude is normalized to unity, we have 

A(t)= (at/n)” exp(n-at) ( 43 ) 

From data obtained from the GSFC experiment on the first Radio Astronomy 
Explorer satellite (RAE-1) we found that this function simulates the rise and 
decay of a typical solar burst observed at a frequency near 0.5 MHz for ns4 and 
a= 1 rain where t measures the time from burst start in minutes. Using Equation 
43 , we proceed with the simulation by introducing computer generated noise and 
spin modulation with characteristics similar to that used for the step signals 
in Figures 4 bo 9- Figure 10 shows the simulation of a burst 30 dB above a 
modulated background. Signal and bacltground for a unity gain antenna with no 
spin are shown in panel A, while panel B shows the observed signal with 
spacecraft spin. Panels C, D, and E show the output of the processing method 
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(with a 12 point prooeasing window) using the data in panel B as input. It is 
clear in panel E how the directions yielded by the processing method switch 
from a value near -60^ to 0* when the burst is present. It is also clear the 
continuously changing amplitude of the simulated burst does not greatly reduce 
the accuracy of the position determination, except near the very beginning of 
the burst where the amplitude change is extremely rapid. 

i|c. Solar Radio Burst .Ql>8ervatlg aa 

The first application of the use spacecraft spin modulation to determine the 
direction of arrival of radio emission was that of Slysh[1967]* By examining 
the phase of observed spin modulation from data maxima and minima of a short 
oblique monopole, Slysh made useful determinations of the arrival direction of 
several solar bursts at 0.2 and 1 MHz. His results were apparently affected by 
limited information on spacecraft aspect. 

The method described in the present paper allows a precise detei’mination to 
be made of the phase by utilizing all the data points in the spin cycle equally . 
In Figure 11, we show an example of a real solar burst measured with the GSFC on 
the IMP-6 spacecraft. The observed modulation is a beat pattern between the 
sampling rate and the spin rate and is similar in duration to the computer 
simulations of Figures 4 through 10. Using the real data at 475 kHz in panel A 
of Figure 11 as input data R(tj^), the processing method outlined in Equations 
30 and 32 yield the results shown in panels B,C, and D respectively. The length 
of the processing Window used was N= 12 points. The derived .source direction 
and source modulation index were then smoothed by a 12 point running average 
before plotting to remove any residual modulation left from a non-uniform 
sampling of phases. Panel D shows the derived source direction. 0® is the 


- 21 - 


direction of the Sun, and positive angles are to the west. Before the onset of 
the burst, the small amount of modulation not readily visible is determined to 
be coming from the direction of the Earth (77“) as seen from the spacecraft. 

At the start of the burst the derived direction of arrival quickly shifted to a 
position about 12® west of the center of the Sun and remained relatively 
constant for the duration of the burst . Note that the derived direction is not 
greatly dependent on the burst signal to background ratio as predicted by 
Equation 37. Panel C is a plot of the derived modulation Index. Before the 
burst, the Earth emission and background have a low level of modulation ( < 

0.1). The burst itself is almost completely modulated (o “ 1) indicating that 
this solar burst at this frequency subtends a small angle from lMP-6 and is 
located very close to spin plane (ecliptic) of the satellite. Panel B is the 
restored burst amplitude using the values of and a (plotted in panels A and B) 
in Equation 32. It is clear that this processing method, applied to actual 
data, yields results similar to those obtained in the simulation shown in Figure 
10 . 

The utilization of these spin modulated techniques on data obtained over a 
wide range of frequencies allowed the first tracking of the course of solar 
radio bursts travelling out from the Sun over distances form 0.1 to 1 AU 
(Fainberg et al,1972|Fainberg,1973;Painberg and Stone, 197^). In a similar way a 
movie was prepared showing the centroid of a solar burst progressing out from 
the Sun along the spiral interplanetary magnetic field in an ecliptic plane 
projection (Fainberg, 1973) . A sequence of frames from this movie is shown in 
Figure 12. 

The use of spin modulation to determine arrival direction has also led to 
the first detection of harmonically related radio emissions (type II) from an 
interplanetary shock propagating out from the Sun at distancnr beyond 5 solar 


radii (Malitson eb al|1973)> 


JJ?.ggoafci,ona Earth BadAo Eaiaaiono 

It is possible to obtain additional information on the location of radio 
sources by triangulation from different directions. This teohniqre is 
particularly effective for satellite observations of radio emissions occurring 
within the magnetosphere of the Earth. Data from the high sensitivity receivers 
of the GSFC experiment on indicate that low level radiation is frequently 

emitted from the Earth at a wide range of frequencies (Stone, 1973) . 

Using the direction of arrival determined form spin modulation, it is 
possible to accurately track the Earth emissions. Figure 13 shows a projection 
of an orbit of IHP-6 (1972) in the ecliptic plane. The directions of arrival of 
radio emission at 130 kHz for selected times are shown as lines from the 
satellite towards the Earth at angles measured from the spin modulation. They 
tend to intersect within a small region on the day side of the Earth’s 
projection. 

By dividing the ^1 day orbit into 3 minute periods, it is possible to form a 
map of over 3 million possible radio arrival direction intersections. Kaiser 
and Stone[1975] have reported the results of such a study using triangulation 
from orbits extending over 500 days of observations in which 2 sources are 
evident. Figure 14 is a summary of their findings at 130 kHz. The source on 
the day side is lower in intensity but is present more often than the second 
(nightside) source which is very sporadic and intense. 
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Mfcefltion Qf Qthep Planetary Emisalona 


By ©xansining the IMP-6 obaepvationa fop peplods when the meaeuped eoupoe 
looationa wep© in the dlpeotion of Jupitep, flpown [197^3 was able to make the 
fipat doteotion of Jupitep padio emlaaion at low frequenoiea. An example of the 
roeaauped dj.peotion of appival of a Jupitep padio bupat ia ahown in Figure 15. 
Prior to the burat which atarted near 0900 on April 25,1972, the meaaured 
direotion of arrival of the radio emiaaion waa weat of the direction of the 
Earth. At the start of the burst, the measured position quickly shifted towards 
the direction of Jupiter, where it stayed for the duration of the burst. In 
this manner. Brown was able to locate, identify, and study a large number of 
Jupiter bursts within the frequency range of 425- 9900 kHz in 1971 and 1972. 

Using similar, direotion finding techniques with the IMP-6 data, Brown [1975] 
was able to make the first detection and measurement of nonthermal radio 
emissions from the planet Saturn. The spectral analysis of these bursts may 
yield valuable information on the planet’s radiation belts and magnetic field 
strength. 

5. SUMMARY AND FUTURE WORK 

The data processing technique described in this paper has allowed a very 
accurate and fast determination to be made of source direotion and character- 
istics using spin modulation of the lMP-6 GSFC radio data. This permitted radio 
techniques to bo used to detect and study processes occuring at long distances 
from the observing spacecraft. In the case of solar emission, it has provided 
the first radio visualization of the curvature of the interplanetary magnetic 
field caused by solar rotation. If several spinning spacecraft separated by 


large dlatanctaa obatnrf the aaite events aouroe looatlon via triangulation from 
extended baaelinea ia poaaiblei independent of modela of emiaaion levela. Suoh 
obaervationa <eg between HELIOS and ISBE 3 ) will alao yield valuable informa- 
tion on groaa magnetic field atruoture and plaama denaity along the trajeotor- 
iea of the partiolea reaponaible for the burat. 

In the ease of emiaaiona from Earth and other planeta, it is clear that the 
radio techniques will permit the exploration via remote aenaing of many 
particle-field interactiona occurring in the planetary magnetospheres. 

Several extensions of these teohniq'ies are possible, For sporadic events 
where the amplitude changes are more rapid, it is possible to derive better 
information on source direction and modulation index by allowing for this. A 
method involving least squares fitting has been derived and utilized on lMP-6 
data and will be de$’oribed in a future report. The least squares technique also 
allows a better way of weighting the data which leads to increased accuracies. 

A very important enhancement of the spin modulation technique is possible if 
the spinning dipole is located out of the spin plane (as it was for Slysh 
CI 967 ]). For that case the spin modulation signal R(t) from a distant source 
for an antenna short compared to wavelength is of the form 

R(t)s [1+ 0 ^ cos(ut-«tQ-e) Og cos2(a)t-utg-4>)3A (44) 

where A is a measure of the source strenth, ♦ is the source direction projected 
into the spin plane, and and Sg are quantities which permit the determination 
of the other angle describing the center of the source direction (elevation 
angle out of the spin plane) . In the case of a uniformly distributed source 
this technique yields both angles describing the centroid of the source, as well 
as the half-width of the source, thus removing the ambiguity illustrated in 
Figure 2. 

A satellite experiment (Knoll et al,1978) based on this latter technique of 
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an oblique apinning dipole Celeotrloally formed ) haa been developed and is 
Inoluded on the Inbernafcional Sun Earth Explorer (ISEE 3, launoiisd in Auguat 
1978) as a joint experiment of the Meudon (France) and OSFC (USA) radio 
astronomy groups. 
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APPENDIX A : ERROR ANALYSIS 


In this appendix exprealona are derived for the errora in the determination 
of the direction of arrival and the modulation index in the presence of an 
unmodulated background. 

1. General Gain Function 

Using the notation as in Section 2, where R(tj^)s observed signal at time 
G(^)s normalized gain in direction A(t^)s signal observed by unity gain 
and B(t^)s background measured at time tj^i we have 

R(t^)= G(^-wt^)A(t^) + B(t^) (45) 

Expanding G in a Fourier series we have 

G(^)s £ bjoosj^+ajslnj^ (46) 

or equivalently 

G(^)= r djcosj(^-ej) (47) 

Assume that the signal ampltude A(t) does not change significantly within the 
processing window of N points and that the mean value of k(t^) over the 
processing window is A. We define the quantities and 

Cj^= 2 R(t^)coskut^ 

Sj^= £ R(tj^)sink(utj|^ (48) 

CqS £ R(t^) 

Since N is a value such that phases of the trigonometric functions appear 
uniform in the sums, we have by standard Fourier techniques the following 
expressions for the mean values 

^k" 2 Abj^cosk(|>g + 2 Aaj^sink*^ 

i At>ij®ink(>g - I Aaj^cosk(>^^ (49) 

Cq= NAbg + NBbg 
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It can be shown that these Equations result from a least squares fitting of 
Equation 45 to the data measured in the processing window. We are interested in 
the variances of Cj^, and Cq, and write these quantities as var var 
and var Cq respectively, We have 

var var C | J djCosj(<}ig-a)t^-ej)A(tj|^)coskut^ ■»-| B(tj^)coakutj^ ] (50) 

Interchanging the i and the j summations 

var Cj^= var [ Z dj J + I (51) 

where 

x^= cosj(<jig-iot^-ej) A(t^) coskut^ 
coskoit^ 

Since x^ and y^ are statistically independent 
2 

var 5 { vaf { var y^ (52) 

Let 0= ut^ and let var^ x^ be the contribution to var for a particular value 
of 0 . Then 

2 2 

var^ x^= cos j(<|)g-0-ej) cos k0 var A(t^) (67) 

Since all 0 (modulo 2tt) are equally probable 

p(0)d0= d0/2ir (54) 

where p is the probability of O being between 0 and 0 + d0. Then 
( 

var x^= j p(0) varg x^ d0 (55) 

Using Equations 57 and 55 for k 0 

var x^= ^ var A for j= 0 

1 

= ^ var A for k 

= ^ [ 1+gcos2k((t)g-ej^) ] var A for j= k/ 0 (56) 

var \ var B 

where var A, var B are the variances of A(tj^) and B(t^) respectively. Using 


Equations 


56 and 57 


vair’ Cj^= I [ dg +2 j dj J oos2k(^g-ej^) ] var A van B 


(57) 


var SjjS i f ^0 j ■ J °o®2k(*g-e|^) ] var A +| var B (58) 

var Cq= 2 C 2dQ + J dj ] var A + var B (59) 

We obtain the standard deviations 
ACj^s /(var C|^) 

ASjjS /(var Sj^) (60) 

ACqS /(var Cq) 

Note that ACj^/Cj^ «. AA/(A/N) « 1 for large signal to background values. This 
applies also to ASj^ and ACq. 


2. Dipole Gain Function 

For a dipole antenna with length short compared to wavelength measuring a 
source with modulation index o we have for the normalized gain function 

G(4>)= (1+acos2(|i)/(1+a) (61) 

From Equations 49 

^2" 5CNAacos24ig)]/( 1+a) 

Sgs 2[NAasin2tg)]/(1+o) 

Cq= (NA+NB)/(Uo) 

From Equations 60 

AC2= /[ 1+ 2 cos4(>g +(1+a)^ AB^/AA^ ] AA /N/[(1+a)/2] 

AS 2 = /[ 1+ 2 cos4tg +(1+a)^ aB^/aA^ ] aA /N/[(1+a)/2] (63) 

ACq= /[ 1+ I a^+ (1+a)^ AB^/AA^ ] aA /N/(1+a) 

To calculate the rms error A^ use 

9 

tan2^g= S 2 /C 2 (64) 

Then 
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1+ I ^ +(1+o)^ AB^/AA^ ] AA/[a/(2N)] 

In a similar way we find the error in a (for large signal to noise case) 

Aa= /[2-3a^/2+aV2] AA/(A/N) (65) 

APPENDIX B - OBSERVED SIGNAL 

In this appendix we derive tne expression for a signal observed by a short 

1 1 

dipole antenna with a gain function G(<}|)= ^ "*■ Section B1 deals with a 

point source located at an angle 6 from the spin axis. Section B2 deals with a 
uniform source. 


B1 


Let ^ be the unit vector pointing toward a point source so that 
^ = (sine^oosi , 3in0„sini)i„, cose.) (67) 

Where 0^ and (ji^ are the polar and azimuthal eingles. Let ^ be the unit vector in 
the spin (azimuth) plane along the dipole axis. Then 

3 . = (cos(ut + 7 t/ 2) , sin(o)t + ir/2),0) (68) 

where ut is the azimuthal angle to the dipole normal in the spin plane. The 
response of the short dipole to the source with strength A(t) is given by 
R(t)= C 1-(a,*^)^ 3 A(t) (69) 

After substitution 

R(t)= [ 1+ aoos2(ii)t-4ig) ] ^ A(t) (70) 

where a is given by 

as sin^0„/(2-sin^0^) 

o o 


( 71 ) 


2. Unilona Conig^i Sounafi 


We consider a source having a uniform intensity over a cone with half angle 
Y and axis located in the y-z plane at an angle D from the dipole axis (z-axls). 
If dfl is a small solid angle within the conical source located at an angle 0 
from the dipole axis, its contribution to the total signal R is 


dRs Asin^odR (72) 

where A is a measure of the source strength. Then 

R= A j sin^edn (73) 

'source 

Initially we place the dipole along the z-axls so that 

Rs A j (1-z^) dn (74) 

' source 


where z is the component of the unit vector towards ds^. To evaluate this 
integral over the source, we first rotate coordinates through an angle 0 about 
the x-axis to a primed set of coordinates such that the negative z-axis is along 
the cone axis. Then 

z= z’cosg-y'sine (75) 


and 

2 2 2 2 2 
z = z’ cos 0+y’ sin 0-2y*z'sinpcosp 

Let O', 'll* be polar cooordinates wiLh respect to the z'-axis. Then 

y'= sinO'cosifr' 

z*= cosO' 

dh= sinOdOd^s sinO'dOd^ 

In the primed coordinates. Equation (74) becomes 

R= A| d*!** I dO* sinO*[i-z'‘^cos^R-y'^ain^e+2y'z'sinpcos0] 

'o 'o 


(76) 


(77) 


(78) 
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( 79 ) 


After substituting Equations (77) we find after integrations that 
R= 2ffA(B-Coo3^0) 

where 

2 1 13 

B= ^ COSY-j oos Y 

1 2 
C= g oosY sin Y 

Equation 79 yields the response of a dipole inclined at an angle P to the cone 
axis. We consider a new xyz coordinate system where the cone axis is the unit 
vector ji, and is inclined at an angle to the z-axis with azimuthal angle in 
the xy plane. We let the dipole unit vector a. be in the xy plane and measure 
the azimuthal angle to its normal from the x-axis. Then 


M = (cos[ut+Tr/23,sinCut+Tr/2]|iO) (80) 

jS. = (sineQOos(ji^,sine^sin(t>Q,cos 0 ^) ( 8 l) 

003 ^ 3 = ( 82 ) 

and finally 

003 ^ 0 = sin^eQ[1-cos2(ut-({ij,)]/2 (83) 

After inserting in Equation 79 we find 

R(t)= Cl+a cos2(a)t-i )]A'(t) (R^l) 

c 

where 

a= sin^e^/( 2 B/C -sin^ 0 „) (85) 

U (J 


and A’(t) is a measure of the source strength. Note that for small cone widths 
Y, Equation 85 reduces to Equation 71. 
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FIGURE CAPTIONS 


1 . 


mUM 2. 


Eisam 3 . 


iifim 4. 


Dipole gain coefficients for the first four terms in the Fourier 
expansion for a range of normalized antenna lengths. For D/X < 0.5» 
only the first two terms are significant. 

Modulation index for a circular conical source with uniform intesity 
as a function of the source width and elevation angle of the source 
center above the spin plane. For example, a modulation index of 0,2 
can result from a narrow source located at an elevation angle near 
55” above ( or below) the spin plane. Wider sources must be located 
at lower elevation angles. 

Two source effects. A) If a second source is present, the direction 
of arrival measured for the combination departs from the direction of 
source 1 by the amount shown. If the second source can be measured 
alone, then a correction can be made analytically. B) If a second 
source is present with modulation index the same as source 1, then 
the combined modulation index measured departs from that of source 1 
by the amount shown. If the second source can be measured alone, then 
a correction can be made analytically. 

Application of the processing method to computer simulated data with 
receiver noise and sampling characteristics similar to those of the 
GSFC experiment on IMP -6. Panel A shows the generated data without 
spin. Panel B shows the generated modulated data with an unmodulated 
background present at all times. At t= 0 min a 10 dB step is switched 
on with a modulation index of 0.75 and with a direction of arrival of 
O”. Panels C, D, and E show results obtained from the processing 
method applied to input data from panel B with a processing window of 
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nquBB 6, 


f-ISUBB 7. 


£IfiUii.E 8. 


miM 9. 


IIQiiM 10. 


12 auooessive points (1 min of data), For the curve marked 
"background removed" the background oharaoteriatios were determined 
by the prooeaaing method uaing 10 min of data prior to the aignal 
atep, and were removed uaing Equation 41. 

Same aa figure 4, but with a aignal amplitude 20 dB above background. 
Background removal ia aeen not to be neceaaary for large aignal to 
background valuea. 

Same aa figure 4, but with a aignal amplitude equal to the 
background. Panel D ahowa the correction provided by the background 
removal techniques (Equations 4l), 

Same as figure 6, but with a processing window 10 times as long. 

Note the improved derived direction of arrival (panel E) and derived 
modulation index ^ panel D) for this example with a processing window 
of 123 points ( 10,5 rain of data). 

Same as figure 7» but with a signal step only 0.1 times background. 
The input signal is nob visible in panel B but, nevertheless, the 
processing method does yield a determination of the arrival direction 
and modulation index. This shows the power of the data processing 
technique for longer processing windows (123 points or 10.5 min of 
data) . 

In this case the background is also modulated with a direction of 
-60° and modulation index 0.25. The signal size is equal to the 
background in amplitude but has an arrival direction of O'® and a 
modulation index of 0.75* The background characteristics are 
measured during a 10 min section of data before the step and are used 
in panels D and E for background removal using Equations 4l. 

Computer generated data simulating the rise and decay of a solar 
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burst. The background is modulated with an arrival direotion of -60** 
and modulation Index 0.75. Panels D and E show the measured 
oharaoteristios to be similar to the burst during the period when the 
burst amplitude Is above the background level. 

FIGURE 11. Actual data from lMP-6 observations of a solar burst measured of 23 
June 1971 are plotted In panel A. Panels B,C,and D, show the results 
of the processing method applied to the data In panel A. The derived 
direction of arrival before the solar burst is determined to be from 
the direction of the Earth (77* )• At the start of the burst the 
direction quickly shifts to a location 12° west of the center of the 
Sun (0°) . 

■ FIGURE 12. Frames from a movie showing the evolution of the centroid of a solar 
burst as determined from the measured arrival direction and 
amplitudes, using lHP-6 data. Each frame is a collection of two 
minutes of data plotted in an ecliptic plane projection with the Sun 
at the center, and the Earth and the IMP-6 satellite at the bottom at 
0°. Circles located at increments of 50 solar radii from the Sun are 
drawn. Emission at different radio frequencies originates at » 

different emission levels in the outer corona of the Sun. The » 

locations of plotted radial line segments (with length proportional J 

to signal amplitude) are determined from the intersections of the 
measured directions of arrival with average spherical coronal 
emission levels. The centroid of a travelling solar radio burst 
caused by energetic particles moving out ft'om the Sun along the 
spiral interplanetary magnetic field is shown evolving in the 
successive 2-minute frames (after Eainberg, 1974) . 

FIGURE 13. Outline of orbit 74 in 1972 of IMP-6 spacecraft; shown are measured 
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arrival direofclons of radio noise ab 130 kHz during seleoted tines 
projected in the eolipbio plane, these direotions intersect in a 
small region on the day side of the Earth. 

FIQUHB 111, Occurence rates of intersections of arrival directions at 130 kHz 
similar to those shown in figure 13. These data summarize 500 days 
of IMP-6 observations (for discussion, see Kaiser and Stone[19753) . 

FIGURE 15. Detection of radio burst coming from Jupiter observed with IMP-6. 0® 
on the angle plot is the direction of the Gun. Positive angles are to 
the west. Prior to 0900, the direction of arrival of the radiation 
quickly shifted towards Jupiter, In this way, the first detection of 
emission from Jupiter at low frequencies was made (Brown, 197^) * At 
this time the beat period between the sampling rate and spin rate was 
much shorter than in figure 11 and is not visible as a spin 
modulation pattern. A processing window of 32 points (2.5 min) was 
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